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Abstract 

An iterative proceedure is proposed to compute the classical gauge field produced in the collision 
of two heavy nuclei at high energy. The leading order is obtained by linearizing the Yang-Mills 
equations in the light-cone gauge, and provides a simple formula for gluon production in nucleus- 
nucleus collisions. At this order fct— factorization breaks down. 



1. Introduction 

Perturbative QCD, based on small coupling calculations, shows that the gluon den- 
sity in the nuclear wave function rises quickly at small x (large energy) and eventually 
saturates to preserve unitarity. The saturation of the gluon distribution is governed by 
a single scale, the so-called saturation scale Qg S> Aqcd, which sets the hard scale for 
pertubative calculations. This idea is formalized in the recently developped theory of the 
color glass condensate fCGO pMHHISI^ . 

In the framework of the CGC, the problem of gluon production reduces to solving the 
classical Yang-Mills equations with a classical current describing the fast moving sources. 
To calculate observables, one has to average over all possible configurations of the sources 
with a given statistical distribution. To probe the saturation regime, proton(deuteron)- 
nucleus collisions have been studied at RHIC since they represent a clean probe of the 
nuclear wave function at high energy. The gluon production cross-section has been calcu- 
lated analytically 7 8 9 10,ll|12j . In this case, the proton is a dilute object, whereas the 
nucleus is saturated. To do so, one linearizes the Yang-Mills equations with respect the 
weak proton field resuming all high density effects in the nucleus, fct— factorization holds 
and allows one to express the cross-section as a convolution in transverse momentum 
space of proton and nucleus gluon distributions. Also, it has been shown that final state 
interactions are absent, therefore no medium is produced. In nucleus-nucleus collisions. 
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where both projectiles are in the saturation regime, numerical studies of gluon produc- 
tion in nucleus- nucleus collisions at high energy, based on the |13ll4j (see also [15]), have 
been performed [16ll7j . but no exact solutions have been found so far. Kovchegov was 
the first to consider analytically the problem, however the formula he proposes is derived 
by conjecturing the absence of final state interactions in addition to other assumptions 
[18] . Later, I. Balitsky calculated the first correction the proton- nucleus by expanding 
the gauge field symmetrically in powers of commutators of Wilson lines [T^] . 

In this work |20| . we propose an analytic method, in the CGC framework, to calculate 
the cross-section for gluon production in nucleus-nucleus collisions. 



2. The classical gauge field in nucleus-nucleus collisions 



The Yang-Mills equations read 

jj^p^u ^ j.^ (2.1) 

where J'^ is a conserved current describing the fast moving nuclei, A and B, moving 
respectively in the -\-z and —z direction. Their source densities p^(a;+,a;) ^ 5{x^)p ^{x) 
and pg{^x~ ,x) ^ S{x~)pg{x) are confined near the light-cone, = {t + z)/\/2 and 
x^ = {t + z)/ \/2 are the light cone variables, and x represents the transverse coordinate. 
As we shall see, the relevant degrees of freedom are the Wilson lines : 

U{x)^Vexp(^ig J dz+ -^p,{z+ ,x) ■ T^, (2.2) 
V{x)^Vcxp(^ig J dz'^Pg{z-,x)-Ty (2.3) 

for nucleus A, and nucleus B respectively. The case where one of the sources is weak 
( corresponding to proton-nucleus collisions) has been solved analytically, however, the 
scattering of two dense systems is still an unsolved problem. In our work we present an 
iterative method to built the solution of the Yang-Mills equation for this problem with 
the aim to capture the main essence of the physics in the first iteration which, as we 
shall see, can bee easily computed. For this purpose, the light-cone gauge = is 
very convenient, and leads to a lot of simplifications. In this gauge the current can be 
constructed exactly, and the field immediately after the collision follows easily from the 
yang-Mills equations. 



The nuclei fields are not determined uniquely, because of the additional gauge degree 
of freedom in axial gauges. One simple configuration is A~ = —-^p^{x~^,x), A^^ = 

= 0, for nucleus A and = — / dy^V^ {y^ ,x) J^p^(j/~,a;) for nucleus B. The 

field just after the collision is determined exactly for this initial condition A^ = UA^^ 
and = A^. The transverse pure gauge field of nucleus B is present all the 

way to i = oo. this would lead to technical complication while computing the field in the 
forward light-cone. To avoid that one can simplify further by removing this pure gauge 
by a gauge rotation involving the gauge link V , i.e., • T ~ V'^{A^^ ■ T)V — j^V^d^V. 
So the produced gauge field near the light-cone gets rotated leading to 
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A'^V{U~l)A'^=ai d+A- = V{dW)A^^=f3o. (2.4) 

Having the exact field produced immediately after the collision of two heavy nuclei one 
can think of solving the equations of motion iteratively in powers of this initial field. Let 
us define an expansion in powers of the initial fields a,, and (3^ : 

QC 

= E ■ (2.5) 

Tl = 

In this case the zero order is simply obtained by gauge rotating the nuclei fields A~ and 

Note that we have generated a + component of the gauge field. Strictly speaking, we 
switch to = — 7p-p^ gauge, which reduces to = in the forward light-cone since 

the source has its support only on the light-cone, — 0. The equations of motion for 
the first order read 

9+(a^4))=0, (2.7) 

□i^,) - d-d,,A^^^ ^ 0, (2.8) 

- d'd^^A^^^ = 0, (2.9) 

which are solved in Fourier space by 



".^A^,(<Z) - -2 - ^ ) aii,) - 2.^/3, (g), 
Vi-^^(g) = -|i4(q). (2.10) 

3. Gluon production 

The spectrum of produced gluon is defined as follows 

i2.f2E^=Y^{\M,\'), (3.11) 

where A is the gluon polarization. The symbol (...) stands for the average over the color 
sources and PbEI]- The amplitude TWa is related to the classical gauge field by the 
reduction formula 

Mx= lim q'A'{q)e\iq), (3.12) 

where f-'xiq) is the polarization vector of the gluon and q its four-momentum. In the axial 
gauge A'^ = 0, only the transverse components of the field contribute, and the sum over 
polarizations states is done with the help of the relation J2^\i^)^*x il) — ■ Note that, 

A 

for on-shell gluons {2q^q^ — = 0), the condition of transversality q^M^^{q) = is 



fulfilled, where = q'^A'^{q). By inserting in (|3.12p the explicit expression of A'^^^^{q) 
given in Eq. (|2.10p we get 

M, = 2 (^4 - ^g^) aliq) - 2z^/3„(g). (3.13) 

This allows us to write the gluon spectrum in the following compact form [20| : 

^^'E^ = X "o(g)P + mq)\')- (3.14) 

Because of the color structure of the fields, it is not possible to write Eq. p.l4p in a 
fct-factorizcd form in the case where the two nuclei are in the saturation regime. In case 
where one of the nuclei is dilute, one can expand Eq. (|3.14[) at first order in the weak 
source (this would be the case for a proton) we recover the well know fct-factorization 
formula for proton- nucleus collisions 111]. 
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